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Abstract 

The neutrino emissivity from two and three flavour quark matter is numerically cal- 
culated and compared with Iwamoto's formula. We find that the calculated emissivity 
is smaller than Iwamoto's result by orders of magnitude when pj(u) +pj(e) —pf(d(s)) is 
comparable with the temperature. We attribute it to the severe restriction imposed by 
momentum conservation on the phase space integral. We obtain an alternate formula 
for the neutrino emissivity which is valid when the quarks and electrons are degenerate 
and p/(u) + Pf(e) — pf(d(s)) is large compared to the temperature. 



It has been conjectured that dense stars may consist of quark matter or quark 
matter core with neutron matter outside [1-5]. Although theoretical understanding of the 
properties of quark matter is not yet available, various quark models have been used to 
calculate the equation of state of the quark matter and determine the properties of quark 
stars [6-9]. Unfortunately, it is found that the properties of quark stars, such as surface 
gravitational redshift z, moment of inertia I, maximum mass M, radius R and pulsar periods 
P, are not significantly different compared to those of neutron stars. Therefore it is difficult 
to distinguish one from the other observationally. 

On the other hand Iwamoto [10,11] has proposed that neutrino emissivity ( e ) could 
play a significant role in distinguishing between quark and neutron stars because it differs 
by orders of magnitude for the two. Particularly e for quark stars is larger by 6-7 orders 
of magnitude than neutron stars which could lead to faster cooling rate for quark stars, 
thus reducing their surface temperature. There are, however, a number of other mechanisms 
[12,13] and modified URCA processes [14] proposed to increase e for neutron matter. 

Iwamoto [10] has derived the formula for e using apparently reasonable approximations 
and this formula has been widely used [8,15-17] to calculate e for two and three flavour quark 
matter. According to his formula e is proportional to baryon density ( Ub ) , strong coupling 
constant ( a c ) and sixth power of temperature (T) for d quark decay. For s quark decay T 
dependence of e is same as that for d decay. Furthermore his results imply that electron and 
quark masses have negligible effect on e and s quark decay ( in case of three flavour quark 
matter ) plays a rather insignificant role. 

In the present paper we want to report an exact numerical calculation of e and a compar- 
ison of our results with the Iwamoto formula. Our results show that the Iwamoto formula 



overestimates e by orders of magnitude when pf(u) +p/(e) —pf(d(s)) is comparable with the 
temperature. For reasonable values of a c and baryon densities, this quantity is much larger 
than the expected temperatures of neutron stars (~ few lOths of MeV) for two flavour quark 
matter, but is comparable with temperature for three flavour quark matter. 
The neutrinos are emitted from the quark matter through reactions 

d — > u + e~ + v e 
u + e~ — > d + u e 
s — > u + e~ + V e 

u + e~ — > s + z/ e (1) 

The equilibrium constitution of the quark matter is determined by its baryon density (%), 
charge neutrality conditions and weak interactions given in Eq.(l). Thus, for two flavour 
quark matter, 

2n u - n d - 3n e = 
n B = (n u + n d )/3 (2) 

and for three flavour quark matter 

Pd = Hu + He{Hv e = Hu e = 0) 
fid = 

2n„ - n d - n s - 3n e = 
n B = (n u + n d + n s )/3. (3) 



The number density of species i is = g.p 3 (i)/ (6n 2 ) with the degeneracy factor ^ being 



two for electron and six for quarks. For electrons /i e = yp 2 (e) + m 2 and for quarks we use 
[18] 

N = [l + S -^(l-^ln( X + V ))] Pf (4) 

where x = pf(q)/m q and i] = Vl + a; 2 , m g being the quark mass. For massless quarks Eq(4) 
reduces to 

N = (1 + (5) 
The neutrino emissivity e for reactions involving d(s) quarks is calculated by using the 
reactions in Eq.(l). In terms of the reaction rates of these equations, we get, 

(Pd(s)-Pu)(Pu-Pe) 



e d ( s ) = A d(s) J d 3 p d(s) d 3 p u d 3 p e d 3 p u - 



E u E d ( s )E e 



x5 4 (p d (s) ~ Pu - Pe ~ Pu)n(p d{s) )[l - n(p u )][l - n(p e )} (6) 

where pi = {E^pi) are the four momenta of the particles, n{pi) = e ^ Ei -^) +1 are the Fermi 
distribution functions and 

24G 2 cos 2 6 C 
24G' 2 sin 2 # c 

A ' = (8 > 
where G is weak coupling constant and 9 C is Cabibbo angle. For degenerate particles, 
(f3pf(i) ^> 1), Iwamoto has evaluated the integrals in Eq.(5) using certain reasonable ap- 
proximations and obtained the simple expressions for e d and e s as given below [11] . 
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td = 315^ ° 0s2 ' dcacP f( d "> p f( u "> p f( e ") T6 



e x = 



457tt 
"840 



G 2 sm 2 6 c a c fi s p f (u)p f (e)T b (9) 



The approximations involved in obtaining these formulas are 

1. neglect of neutrino momentum in momentum conservation, 

2. replacing the matrix elements by some angle averaged value, 
and 

3. decoupling momentum and angle integrals. 

The expressions for neutrino emissivity as obtained by Iwamoto have been used widely. The 
temperature dependence of emissivity as obtained by Iwamoto has a physical explanation. 
Each degenerate fermion gives one power of T from the phase space integral ( d 3 pi — > 
Pf(i) 2 dEi oc T ). Thus one gets T 3 from quarks and electrons. Phase space integral for the 
neutrino gives d 3 p u oc (El)dE u oc T 3 . Energy conserving 8— function gives one T" 1 which is 
cancelled by one E v oc T factor coming from matrix element. So finally one gets e oc T 6 . 
This argument, however, ignores the fact that Apa (Ap s ) = Pf(u) +Pf(e) — pf(d) 
which is related to the angle between f>d , p u and p e could be small and comparable to T. 
We shall demonstrate below that precisely in this region that the Iwamoto formula fails. 

Before discussing the causes of the failure of Iwamoto formula, let us first compare our 
results with the Iwamoto formula and try to find out the specific cases where the deviation 
is more pronounced. In Figs. 1-3 we have plotted e vs T for 2-flavour d decay, 3-flavour d 
decay and s decay respectively. For 2-flavour d decay our results (e<z) are in good agreement 
with the emissivity calculated using Iwamoto formula (^d(s)i)- In Fig.l curves (a) and (b) 
are and eai respectively, for a c = 0.1 and n# = 0.4. (c) and (d) corresponds to the 
same but for a c = 0.1 and rib= 1.4. It is evident from the figure that agreement of Iwamoto 



results with our calculation is better for higher densities and lower temperatures. Also td is 
consistently smaller than e^i, the Iwamoto result, in the range of temperatures considered. 
Corresponding fermi momenta of quarks and electrons are given in Table 1. It is to be noted 
that all the momenta are much larger than the temperature. 

Fig. 2 shows the td for 3-flavour quark matter. It shows that €di is 2 -3 orders of magnitude 
higher compared to our results. Here , contrary to the two flavour case, the difference becomes 
more pronounced at higher densities. Fig.3 shows the variation of e s with temperature. Here 
again it is clear that e s is quite different from e s j but this difference is less compared to 
that between and e^i- For all the cases the difference between our results and those 
using Iwamoto formula increases at higher temperatures. The fermi momenta of quarks and 
electron for 3- flavour case are given in Table 2. The study of all the figures and tables above 
reveals that the cases where Iwamoto formula agrees reasonably well with our results, Ap d 
( or Ap s ) is much larger than the temperature. On the other hand when this difference is 
smaller or comparable with the temperature, the Iwamoto formula overestimates the exact 
result by order of magnitude. In addition to these, Fig.3 also shows that our results are 
about a factor of 2.5 lower than the Iwamoto results even at lower temperatures. We have 
found that this difference comes from the approximation involved in the calculation of matrix 
element. 

Furthermore Table 2. shows that for three flavour case electron chemical potential ( 
which is same as p/(e) for massless electrons ) becomes small ( < l.MeV ) for some values 
of a c , riB and m s . In these cases, electrons are nolonger degenerate. Clearly, for such cases 
the Iwamoto formula is not applicable. This point is missed in earlier calculations. 

Our results have profound implications on neutrino emissivity and quark star cooling rates 



because all the earlier calculations have used the Iwamoto formula and predicted large quark 
star cooling rates in comparison with the neutron star cooling rates for temperatures less 
than 1 MeV. Our results show that, particularly for 3-flavour quark matter, the calculated 
emissivity is at least two orders of magnitude smaller than the one given by Iwamoto formula 
and therefore, the three-flavour quark star cooling rates are that much smaller. Hence it is 
necessary to understand why Iwamoto formula fails. 

To investigate the failure of the Iwamoto formula, we consider the integral 



r d p d d p u d p e d gffi ^ _ pu _ pe - p v ) n (p d{s) )[l - n(p u )\[l - n(p e )\. (10) 
J e d e u e e 



Here, we have replaced the neutrino emission rate by unity and therefore I is essentially the 
phase space integral. Following the reasoning of Iwamoto, this integral should be proportional 
to T 5 . Choosing the coordinate axes such that p d is along z-axis and p u is in x-z plane and 
using the 3-momentum 5— function to perform electron and u-quark angle integrations, we 
get 



j = 8% 2 I PddpdVldpuVldped^Pu \j l ~ x l 



e d e u e e p e p u (^Jl - x\{p d - p v x v ) + p v x v ^j\. - x 2 v cos (j) u 

S(e d - e u - e e - e v )n(p d )[l - n(p u )\[l - n{p e )\. (11) 



where x u = cos 6 U and x u = cos 9 U is determined by solving 

PuXu = Pd-PuXu - \p 2 e - pl(l - xl) - pl(l - xl) 



- 2p u p^(l - xl)(l - xl) cos^J 1 / 2 . (12) 

The integral in eq(ll) above is restricted to the momenta \pi — Pf{i)\ few times T due 
to Fermi distribution functions and the energy 5— function. Now, if we neglect the neutrino 



momentum in the 5— functions, we get, x u = (p d + p\ — pl)/2pdPu and the factor in the 
square brackets of eq(ll) becomes l/p d p u p e . 
Two points should be noted at this stage. 

1. Generally, x u is close to unity, so that 1 — x\ is small. But, if Ap d is of the order of 
T, \Jl — %uPd can be comparable with T and p v and therefore p v cannot be neglected 
in the momentum 5— functions. Particularly, the denominator in the square bracket 
of eq(ll) cannot be approximated by p e PuPd\j^- — Thus, if Pd\J^ — %\ < Pv, one 
would get a power of T from the denominator and I will not be proportional to T 5 . 

2. Secondly, the momenta may differ from the corresponding Fermi momenta by few times 
T in the integral. When Ap d ~ T, there are regions in PdPuPe— space where x u > 1 
and the rest of the integrand is not small. Clearly, these regions must be excluded 
from the integration as these values of x u are unphysical. If one does not put this 
restriction, as is done when one factorises angle and momentum integrals, the phase 
space integral will be overestimated (and wrong) when Ap d ~ T. 

The above discussion clearly shows why the integral in eq(ll) should not be proportional 
to T 5 when Ap d ~ T. In order to demonstrate this point, we have calculated the integral 
in eq(ll) numerically and compared with the approximation where the neutrino momenta 
are neglected and the restriction imposed by x u condition is not imposed. The calculation 
is done for a c = 0.1 and for two-flavour case. The results are shown in Fig(4). In this figure 
we also show the result for a case where the electron mass is taken to be 25 MeV. This is of 
course unphysical, but by adjusting the electron mass we can reduce Ap d . The figure clearly 
shows that the approximate value of / is proportional to T 5 where as the exact integral is 



smaller than the approximate one at large T. Further more, for 25 MeV electron mass, the 
departure from T 5 sets in at smaller value of the temperature. This clearly shows that the 
departure is dependent on the value of Ap d . Here we would like to mention that for some 
values of a c and m s , p/(e) is small and is of the order of T. This implies that electrons are 
no longer degenerate and deviation from the Iwamoto result is most pronounced. 

In eq(ll) we have dropped the matrix element of the weak interaction in the emissivity 
calculation (eq(6)). So, the discussion of preceeding paragraphs apply to the emissivity 
calculation as well. Therefore it is now clear why the Iwamoto formula fails when Ap d (or 
Ap s in case of weak interactions involving strange quarks) is close to the temperature of the 
quark matter. We would like to note that the failure of T 6 dependence of the emissivity 
essentially arises from invalid approximations in the phase space integration. Since similar 
arguments are used to obtain the neutrino emissivity of neutron matter, it is possible that 
the emissivity calculated for neutron matter may also be overestimated. We are investigating 
this point. 

Since the departure from the Iwamoto formula arises from the fact that T / Ap d (or T/Ap s 
) is not small, it may be possible to fit the numerically calculated e with a function of the 
form eif(x), where x = T / Ap d (T '/ ' Ap s for strange sector ). The function f(x) should be 
such that for small values of x it should approach unity. Choosing fix) = -7— , 1 , 2 - — 3, 

_L . ~x~ CLX ~\~ OX ~\~ CX 

we have fitted the calculated e for a number of values of a c and m s and obtained the 
values of a, b and c. The quality of fit is shown in Fig. 5. The values of a, b and c are 
—2.0, 110. and 30. respectively. Here we would like to point out that since, for s- decay, as 
mentioned above, there is a difference of a factor of 2.5 in Iwamoto and our results even at 
lower temperatures, the data points for s- decay, in Fig. 5., have been scaled accordingly. 



To summarise, we have demonstrated that the Iwamoto formula of neutrino emissivity 
fails when T/Ap d (or T/ Ap s ) is not small. The formula fails because the neglect of neutrino 
momentum and factorisation of angle and momentum integrals is not valid. We propose an 
alternate formula which is obtained by fitting the numerically calculated e. 
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FIGURE CAPTIONS 
FIG.l. Two flavour <i-decay for a c = 0.1; (a) Iwamoto results for n B = 0.4/m~ 3 , (b) Our 
results for n B = 0Afm~ 3 (Ap d = 6.78),(c) Iwamoto results for n B = lAfm~ 3 , (b) Our 
results for n B = lAfm~ 3 (Ap d = 10.29). 

FIG. 2. Three flavour (i-decay for a c = 0.1 and s quark mass is 150 MeV; (a) Iwamoto 

results for n B = lAfm~ 3 , (b) Our results for n B = 1.4/m~ 3 (Ap^ = 0.067),(c) Iwamoto 

results for n B = 0.4/m~ 3 , (b) Our results for n B = 0.4/m -3 (Ap d = 0.39). 

FIG. 3. Three flavour s-decay for a c = 0.1 and s quark mass is 150 MeV; (a) Iwamoto 

results for n B = lAfrn" 3 , (b) Our results for n B = lAfm" 3 (Ap s = 1.613),(c) Iwamoto 

results for n B = 0Afm~ 3 , (b) Our results for n B = 0.4/m" 3 (Ap d = 9.719). 

FIG. 4. Two flavour phase space integrals for a c = 0.1 (a) Without restriction on cos# u for 

both electron mass m e =0.0 and 25 MeV, (b) Exact integral for m e =0.0, (c) Exact integral 

for m e =25 MeV. 

FIG. 5. ^ is plotted aganist x where x = - . The fitted function is f(x) = 

1. + ax + bx 2 + cx 3 where a = -2.0, b = 110. and c = 30. 



Table 1. Baryon number density ub, Fermi 
electron p/(e) for different a c , where Ap d = 



momenta of u-quark pf(u) 
Pf( u ) + Pf( e ) -Pf( d )- 



, d-quark Pf(d), and 





n B 


p f (u) 


Pf(d) 


P/(e) 


Ap d 




(fm- 3 ) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 




0.60 


357.80 


449.20 


99.15 


7.75 


0.1 


f.00 


424.22 


532.52 


117.56 


9.20 




1.40 


474.57 


595.79 


131.51 


10.29 




0.60 


357.71 


449.25 


95.43 


3.89 


.05 


f.00 


424.11 


532.65 


113.15 


4.61 




1.40 


474.45 


595.87 


126.57 


5.15 



Table 2. Baryon number density Ub, Fermi momenta of u-quark ]?/(«), d-quarkp/(d), s-quark 
Pf(s) and electron p/(e) for different m s and different a c , where Ap d — Pf(u)+pf(e) —pf(d) 
and Ap s = p f (u) + p f (e) - p f (s) 



m s 


a c 


n B 


p f (u) 


p f (d) 


Pf(s) 


P/(e) 


A Pd 


A Ps 


(MeV) 




(fm' 3 ) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 






0.60 


356.99 


360.06 


353.86 


3.33 


0.26 


6.46 




0.1 


1.00 


423.26 


424.81 


421.69 


1.69 


0.14 


3.26 


150.0 




1.40 


473.49 


474.27 


472.72 


0.84 


0.06 


1.61 






0.60 


356.99 


365.89 


347.62 


9.28 


0.38 


18.65 




0.05 


1.00 


423.26 


430.29 


415.98 


7.33 


0.30 


14.61 






1.40 


473.49 


479.49 


467.34 


6.25 


0.25 


12.40 






0.60 


356.99 


365.93 


347.58 


9.70 


0.76 


19.11 


200.0 


0.1 


1.00 


423.26 


429.10 


417.25 


6.34 


0.50 


12.35 






1.40 


473.49 


477.66 


469.25 


4.52 


0.35 


8.76 






0.60 


356.99 


374.26 


337.85 


18.00 


0.73 


37.14 




0.05 


1.00 


423.26 


437.14 


408.39 


14.47 


0.59 


29.34 






1.40 


473.49 


485.45 


460.90 


12.46 


0.50 


25.05 



